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Using the three-site Hubbard model and Floquet theorem, we investigate the 
dynamical behaviors of two electrons confined in a line-shape three quantum dot 
molecule driven by an ac electric field. Since the Hamiltonian contains no spin- 
flip terms, the 6 dimension singlet state and 9 dimension triplet state sub-spaces 
are decoupled and can be discussed respectively. In particular the 9 dimension 
triplet state sub-spaces can also be divided into 3 three-dimensional state sub-space 
which are fully decoupled. The analysis shows that the Hamiltonian in each three- 
dimensional triplet state sub-space as well as the singlet state sub-space in the no 
double-occupancy case has the same form which is similar to that of the driven two 
electron in two quantum dot molecule. By solving the time-dependent Schddinger 
equation, we investigate the dynamical properties in the singlet state sub-space, 
and find that the two electrons can maintain its initial localized state when driving 
by an appropriately ac-field. In special, we find the electron interaction enhances 
the dynamical localization effect. Using both perturbation analytic and numerical 
approach to solve the Floquet function leads to a detail understanding of this effect. 



PACS: 72.20.Ht; 78.66.-w; 33.80.Bc 



I. INTRODUCTION 

Quantum state engineering via optical or electrical manipulation over the coherent dy- 
namical of suitable quantum mechanical systems is a subject of great current interest because 
of a growing number of possible experimental applications. The development of laser and 
masers also open the doorway for creation of novel effects in nonlinear quantum systems 
which interact with strong electromagnetic field^. In particular, recent experimental sue- 
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cesses in detecting Rabi oscillations in quantum dot(QD) systems driven by ac- field has also 
spurred interest in the use of driven ac-field to coherently manipulate the time development 
of electronic states^. Achieving an understanding of the dynamics of the QD system is 
extremely desirable, as the ability to rapidly control the localization of electrons suggests 
possible application to quantum computation and quantum information processing, in which 
the coherent manipulation of entangled quantum states is essential component. An exciting 
possibility is to make use the phenomenon of coherent destruction of tunnelling (CDT), 
in which the tunnelling dynamics of a quantum system become suppressed at certain pa- 
rameters of the field. Tuning the driving field thus provides a simple mechanism to localize 
or move charge within the quantum dot on a rapid time-scale by destroying or restoring the 
tunnelling between regions of the devices, so allowing the ac-field to be used as "electron 
tweezers"—. 

In this paper, we address the dynamical localization of two interacting electrons con- 
fined in three QD molecule with linear arrangement driven by an ac electric field. We 
use an three-site model of Hubbard-type to describe the dynamical system, which gives 
a considerable computational advantage over standard numerical approaches, and also al- 
lows us to easily include the important effects of the electron correlations produced by the 
Coulomb interaction. The two electron distribution states are taken as basis vectors to 
write the Hamiltonian matrix. We show that the system can be analyzed in singlet state 
and triplet state sub-space respectively, and the triplet state sub-space can also be divided 
into 3 three-dimensional syate sub-space. By solving the time-dependent Schodinger equa- 
tion, we investigate the dynamical properties in the singlet state and triplet state subspace, 
and explain these finding by Floquet theorem. These exposition are useful to understand 
how the ac electronic field can affect the charge distribution inside a QD molecule, and how 
they can be used for quantum control. 

II. MODEL AND METHODS 

We consider a highly simplified model in which each quantum dot is replaced by a single 
site. The electron can tunnel between in the sites, and importantly, we include interactions 
by means of a Hubbard [/-term. We do not take the energy level detuning between the 
three quantum dots into account. Assuming that the AC field frequency is much lower that 
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single-particle level spacing, the higher-lying single-particle states are also ignored. The 
Hamiltonian defining the system read: 

H(t) = E e*k(t)dl k d ak + £ W(dl k d ak+1 + hc) + 

(er=U)(fc=l,2,3) (<T=U)(As=l,2) ^ 

E U^n^n^ + n ifc n Tfc ) + E ^2(n CTfe n CT ' fjfc+1) ). 

(fc=l,2,3) ( CTl( j')(fc=l,2) 

Where cr is the electron spin, and f and j are used to indicated the up-spin and down- 
spin; k = 1,2, 3 is the serial number of the QD; d} ak jd ak are the creation/annihilation for 
an electron of spin a on dot k; The quantity W denotes the hopping between adjacent dot; 
Ui is the standard Hubbard-f/ terms, giving the energy cost for the double-occupation of 
a dot; Ui presents the Coulomb repulsion between electrons occupying neighboring dots; 
n ak — d ak d a k is the total charge occupation of dot k. For convenience, we consider the three 
QD molecule is in a standing wave field, and the centra QD is just in the wave node. We 
also let the QD molecule align with propagation direction of the standing wave field. So the 
energy level can be parameterized as: e a k(t) = ±(fc — 2)Vcos(ut), where u is the frequency 
of driving field, and V is the amplitude of energy level. Since V is proportionate to the 
amplitude of the ac-field, we still take V as the amplitude of the ac-field conveniently. 

There are 15 kinds of two electron distribution in the three QD molecule which are shown 
in Table HI where A,B,C indicate the site of the three QDs. We replace the states |5), |6), 
|9>, |10), |13), |14> by |5'> = ^(|5> + |6», |6'> = ^(|5> - |6», |9'> = ^(|9> + |10», 
1 10') = 75(|9) - |10», 1 13') = ^=(|13> + |14)), \U') = ^(|13) - |14», and take them as 
basis vectors. A major advantage of these is that they manifest the two-electron spin and 
space distribution clearly. We can write the Hamiltonian(Eq^) in the space spanned by 
the basis vectors. Therefore the Hamiltonian(Eq^) can be described by 15 x 15 matrix. 
The Hamiltonian(Eq^) contains no spin-flip terms since measurements on semiconductor 
QD show that the spin-flip relaxation time is typically extremely long, and so the singlet 
state and triplet state sub-space are completely decoupled. Thus we can investigate the 
dynamical behaviors of two electrons in two sub-space respectively. Therefore, if the initial 
state posses a definite party this will be retained throughout its time evolution, and we only 
need to include states of the same party in the basis. 

The Hamiltonian(Eq^) also contains no spin-orbital coupled terms, and so the wave 
function can be written as the direct product of spin and orbital wave function. A many- 
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particle basis can then be constructed by taking Slater determinants of single particle states 
defined on the QD. The basis vectors (|4), |12), |8)) constitute the up-spin sub-space, and 
the basis vectors (|7), |15), |11)) constitute the down-spin sub-space. Their spin wave 
function are symmetric under particle exchange. The orbital wave function of basis vector 
(|6'), |14), 1 10 ) ) are antisymmetric under particle exchange, so their spin wave function 
are symmetric. Therefore, the basis vectors (|4), |12), |8), |7), |15), |11), |6'), |14'), |10')) 
constitute 9-dimension spin-triplet state sub-space. The calculation also shows that the 3 
three-dimensional spin symmetric sub-space are completely decoupled, and their Hamilto- 
nian spanned by three spin-symmetric basis vectors hold the same form which is 



( U 2 - V cos(cut) W \ 



H(t) 



(2) 



WOW 
\ W U 2 + Vcos(ut) J 

So the dynamical behaviors of the two electrons are identical in the 3 three-dimensional 
sub-space. For simple, we only investigate the dynamical characters of two electrons in 
up-spin sub-space (|4), |12), |8)). 

In order to obtain the dynamical characters of the three-level system, it is necessary 
to examine the system of two electrons confined by two QD molecule driven by ac-field, 
which has been fully investigated in Ref . 4 i 5 i 6 i 9 . In this system, the basis vectors 
(|1, 1) — |1, l))/v^2, |1, 1)) constitute the triplet stste sub-space, and the basis vectors (|2,0), 
(|1,1) + |l,l))/\/2, 1 , 2) ) constitute the singlet state sub-space. The two sub-space are 
also decoupled. In the triplet stste sub-space, the two-particle basis vectors (|1, 1), (|I, 1) — 
|1, 1))/v2, 1 1 , 1) ) are the eigenvectors of its Hamiltonian and constitute the trivial triplet 
state sub-space in which the electron number on each QD is invariably one, and the time- 
dependent term does not influence this characteristics. In the singlet state sub-space, the 
Hamiltonian is similar to that of Formula (EqJ2J). The relationships of the parameters and 
the basis vectors between the two system are: U 2 — > U\ — U 2 , W — > \^2W, |4) — > |2,0), 
|8)^|0,2), |12)^(|l,l) + |l,l))/v / 2. 
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III. RESULTS 

A. The two-electron dynamical properties in the triplet state subspace 

We now investigate the dynamical properties of the three-dimensional triplet subspace 
and take the subspace (|4), |12), |8)) for example. The time periodicity of Hamiltonian(Eq|2J) 
enables us to describes the dynamics within the Floquet formalism. In addition, since 
the Hamiltonian is invariant under the combined dynamic parity operation z — > —z; t — > 
t + Tr/u, each Floquet state is either odd or even. Quasi-energies of different party may 
cross, otherwise an avoided crossing may occurs^. We calculate the quasi-energies and P m in 
as function of V with the parameters U2 = nu. We used two method to obtain quasi-energies: 
(1) the numerical method^ which is to diagonalize the unitary time-evolution operator for 
one period of the driving field U(t+T, t) . It may be easily shown that the eigenvectors of this 
operator are equal to the Floquet states, and its eigenvalues are related to the quasi-energies 
via Xj = exp(-iejT). This method is particularly well-suited to our approach, as U(T, 0) 
can be obtained by integrating the unit matrix in time over one period of the field using the 
Runge-Kutta method. (2) the perturbation approach^. In this method, the quasi-energies 
are obtained by first solving the Floquet equation^ in the absence of tunnelling terms, and 
then performing perturbation theory with the tunnelling terms as the "perturbation". The 
perturbation solution are: E\p = ±\/2J n (V J 'oS)W , £3 = 0. where, J n is the bessel function 
of the n order. We now define the probability P m in- In this singlet state sub-space, the 
wave function can be written as a superposition of basis vectors: \ip(t)) = J2C a (t)\a), \a) 

a 

is the basis vector and a = 4, 12, 8. Substituting this expansion into the time- dependent 
Schrodinger equation yields a first-order differential equation for the expansion coefficients 

c a (h = iy. 



.d 



\C 4 J 



H(t) 



( cA 



(3) 



When the initial condition C a (0) are given, we use the Runge-kutta method to solve the 
differential equation (EqEJ) and obtain the expansion coefficients C a (t). When 6*4(0) = 1, 
£12(0) = 0, C 8 (0) = 0, We term the minimum value of |C 4 (t)| 2 of evolving 25 driving period 
LP min . When C 4 (0) = 0, C 12 (0) = 1, C 8 (0) = 0, We term the minimum value of |C 12 (t)| 2 
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of evolving 25 driving period CP min . LP min and CP min are used to quantify the degree 
of quantum tunnelling. If LP min = 1 and CP m i n = 1, it means that the two electrons 
can maintain its initial localization state in a short time span, and the tunnelling between 
different electron distribution is suppressed completely. If P m in = (or CP m i n = 0), it 
means that the two electrons can not maintain its initial state. 

We calculate the quasi-energies s and the quantity LP min and CP min in FigHJas a function 
of V with the parameters U2 — 8, W — 1, u — 1. The calculation show that the excellent 
agreement between the quasi-energies exact solution and perturbation solution for strong 
and moderate fields. For weak fields, however, the driving terms do not dominate the 
tunnelling terms and the perturbation theory breaks down. Figure d also shows that the 
dynamical properties can be divided into two regime: the strong field regime and the weak 
field regime. In the strong field regime, the peaks in LP min and CP m i n is found by locating 
the roots of J%(V/uj). So the initial states |4), 1 12) and |8) can be maintained. Hence the 
tunnelling between different electric distributions is depressed. In the weak field, CP m i n 
decays smoothly to zero, and a series of peak of LP m i n is found at the quasi-energy crossing. 
So the initial state 1 12) fail to hold on when increasing the driving ac field. However the 
initial state |4) and |8) keep its initial state even though the strong interaction between 
them. 

To interpret the dynamical behaviors for the strong field regime, we seek the perturbation 
approach^. We treat the tunnelling terms in Hamiltonian (Eqj2j as a perturbation. Solving 
the Floquet equation, the degeneracy Floquet states of the zeroth-order approximation are 
given by (U2 = ncu): 

= I exp[— iU^t + %¥- sinatf], 0, 0} 
< MO) = |0,1,0) (4) 
1^3 (^) ) = |0, 0, exp[— i[/ 2 t + i^- sincui]) 
With the degeneracy Floquet states, we obtained the perturbing operator P^ = 

((ui(t)\H t \uj(t))) , where ((■ ■ •)) denotes the inner product in the extended Hilbert spaced 

00 

By using the identity: exp[— iftsmuit] = Jm{@) exp[— imut] to rewrite the form of 

m=— 00 

\ui(t)), the matrix element of P can be found to be P\2 = P23 = W^^n(^)- When the ratio of 
the field strength to the frequency is a root of Bessel function J n , the transition between two 
Floquet states |4) and 1 12) as well as 1 12) and |8) is forbidden. Therefore, the tunnelling be- 
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tween different electrical distribution is depressed intensively, and the phenomenon of CDT 
occurs. 

Let us analyze the dynamical properties in the weak field regime. In order to show it 
clearly, we show in Fig|2^a) the magnified view of Fig^a) for the weak field regime. The 
FigEfa) shows that the dynamical localization occurs at the crossing of quasi-energies e 2 
and e 2 . In the absent of an external driving field, the eigenvalues and eigenvectors are 

' \y x ) = \l,-a/(V2W),l>, E 1 = b= {-yful + 16W 2 + U 2 )/2 
< \<p 2 ) = I -1,0,1 >, E 2 = U 2 (5) 

|^ 3 > = |l,-6/(V2W0,l>, E 3 = a = (yJu2 + 16W* + U 2 )/2 
When driving by weak field, the two-electron Floquet states are approximated by their 
eigen-states. We note U 2 ^> W, and the eigen-state | <^ 3 > has small component of two- 
electron state 1 12), so the eigen-state \(p 2 > and \(ps > are similar to the excited state 
(asymmetry) and the ground state (symmetry) of a electron in two QD. Thus it expects 
the initial two-electron state |4) and |8) remains forever at the crossing of quasi-energies e 2 
and £ 3 4 i 5 i 9 . The rate V/oj at the first crossing is about 1.2, a root of the zero-order Bessel 
function, suggesting that this kind of dynamical properties can be approximates by a two- 
level model. When increasing the interaction U 2 , the component of state 1 12) in |</? 3 ) will 
decrease. So the Coulomb repulsion may help to maintain the initial states |4) and |8). We 
present the calculation in Fig|2fb) with U 2 = 16, and the other parameters is the same as 
that of Fig|2fa). The calculation in Fig|2fb) shows that a series peaks of LP min emerge, and 
its height increase obviously with the increasing of U 2 . 

B. The two-electron dynamical properties in the singlet state subspace 

We now investigate the dynamical behaviors for the singlet state sub-space. The orbital 
wave function of the basis vectors (|1), |2), |3)) are symmetric under the electron exchange, so 
their spin wave function are antisymmetric. The spin wave function of the basis vectors (|9), 
1 5 ), 1 13 )) are antisymmetric under the electron exchange, so their orbital wave function are 
symmetric. Therefore, the basis vectors (|9'), |2), |5'), |3), |13'), |1)) constitute the singlet 
sub-space. The HamiltonianQ spanned by the six spin-antisymmetric basis vectors can be 
written as 6 x 6 matrix, which is 



s 



H(t) 



U 2 + Vcos{ut) 

V2W Ui 
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V2W U 2 + Vcos(cut) 



V2W 





w 
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V2W 



U 1 + 2Vcos(ut) 







Vzw 










Ui - 2Vcos(ujt) J 



(6) 



Due to the diversity of the size and couple between the QDs, the dynamical behaviors of 
the system show multiplicity. We begin our investigation by first considering the simplest 
case, that of the Hubbard L^-term to be infinitely large - that is, we work in the sub-space 
of states with no double occupation. Our Hilbert space is thus three-dimensional spanned 
by the basis vectors (|5'), |13'), |9 )). Therefore, the Hamiltonian can be written as: 

' U 2 -Vcos(ujt) W N 

H(t) W W (7) 

v W U 2 + V cos(W) j 

This form of Hamiltonian(EqEJ) is the same as that of Hamiltonian(Eq|2I). Thus the 
six-level system can be simplified to the three-level system, and the dynamical characters in 
this special case is also similar to two electrons in two QD molecule. 

We now take the most general case, and consider the U\ and U 2 to be a finite value 
which means that the three double-occupation are no longer energetically exclude from the 
dynamics, and accordingly we must take the full six-dimensional basis set. we calculate the 
quasi-energies and P m in as function of V with the parameters U\ = rtiu and Ui = n 2 u>. We 
used two method to obtain quasi-energies: (1) the numerical method; (2) the perturbation 
approach, which are similar to the calculation of Figure ^ The perturbation solution of 
quasi-energies are 



£l,2 

£3,4 
£5 



±V2J ni . n2 W 







6J 2 W 

<Jt -'ni— ri2 ' r 



(8) 



where, J ni _„ 2 and J n2 are the bessel function of the ri\ — n 2 and n 2 order respectively. 
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We now define the probability P m in in the singlet sub-space. The wave function can be 
written as a superposition of basis vectors: \i/)(t)} = C a (t)\a), \a) is the basis vector and 

a 

a — 9 , • • •, 1. Substituting this expansion into the time- dependent Schrodinger equation 
yields a first-order differential equation for the expansion coefficients C a (Ti = 1): 



(c,\ 



.d_ 

'dt 



H(t) 



(9) 



V C x j \ ClJ 
When the initial condition C a (0) are given, we solve the differential equation(EqlHI) by 
using the Runge-kutta method to obtain the expansion coefficients C a (t). We term the 
minimum value of |C a (t)| 2 of evolving 25 driving period P m i n , and use it to quantify the 
degree to which the ac-field bring about localization. In special, we term the minimum value 
of |C*i(t)| 2 of evolving 25 driving period AP min when Ci(0) = 1, C Q (0) = a ^ 1. We also 
term the minimum value of | C2 (^) | 2 of evolving 25 driving period BP m i n with the initial 
condition 6*2(0) = 1, C a (0) = a ^ 2. If AP min = 1, it means that the two electrons can 
maintain its initial localization state (the left QD A) in a short time span, and the tunnelling 
between different electron distribution is suppressed completely, so the phenomenon occurs. 
If AP m i n = 0, it means that the two electrons can not maintain its initial localization state 
(the left QD). 

In order to show the dynamical properties in this general case, it is necessary to show 
the case of non-interacting and double-occupancy permitted firstly, and then investigate 
how the Coulomb interaction affects its dynamical behaviors. The three sites model in 
this non-interaction case is similar to that non-interaction electrons driven by ac-field in a 
superlattice^, with the quasi-energy crossing corresponding to "miniband collapse". We 
term the phenomenon as "quasi-energy collapse". The electron tunnelling between QD is 
considerably suppressed with parameters at the roots of Jo(V/u) (Jo is the Bessel function 
of first kind). The calculation in Figure Efa) verifies our assertion. We show in Figure Hit a,) 
the quasi-energy spectrum and AP min obtained by sweeping over V. The dot and slimline 
present the exact solution and the perturbation solution of quasi-energies respectively, and 
the thick line presents AP min . Figure Efa) demonstrates the excellent agreement between 



10 



the exact solution and perturbation solution of the system quasi-energies for strong and 
moderate fields. For weak fields, however, the driving terms do not dominate the tunnelling 
terms and the perturbation theory breaks down. Figure a) also shows that the peaks in 
AP m in by locating at the point of quasi-energies crossing. Hence the phenomenon of CDT 
occur at the point of quasi- crossing. 

When the electron interactions involving, the difference from the non-interaction case 
arise. Using the same method as that of Figure a), we calculate in Figure Efb) the quasi- 
energy and AP min as a function of V with U\ = 12 and U2 = 4. It shows that the phenomena 
which is similar to " miniband collapse" occurs in the strong-field regime of the quasi-energy, 
even though not all the quasi-energy levels cross at one point. At the weak-field regime, the 
feature of quasi-energy (dotted line) is different obviously from that of the non-interaction 
case, and a series of quasi-energy crossing appears. The calculation in Figure Efb) also shows 
that y\P m i n V can be divided apparently into two regime: the weak-field regime and the 
strong-field regime. At the strong-field regime, the character of AP mm — V is similar to that 
of the non-interaction case (FigO^a)). At the weak field regime, it is interesting to see that 
some sharp peaks of P m in occur at the quasi-energy crossing. The results indicate that two 
electrons localized in left QD initially can maintain its localization state (a 25 driving period) 
when driving by a weak-field, even though there is a strong Coulomb interaction. Therefore, 
the Coulomb interactions have little effect on strong-field regime, and they enhance the 
phenomena of CDT in the weak-field regime. This quantum phenomena can also occur in 
two electrons in two QD molecule, and have been found and analyzed in Ref^. 

We now analyze the system dynamical properties when driving by weak field. In order to 
show it clearly, we give in Fig|3fc) the magnified view of FigCJfb) for the weak field regime. 
FigEHb) demonstrates that the position of the peak in LP min to be found by locating at 
the crossing of quasi-energies E\ and Ei- When driving by weak field, the Floquet states 
can be approximated by the corresponding unperturbed eigen-states. We diagonalize the 
Hamiltonian (EqEJ) when V = 0, and obtain ^ = (-0.5774, 0.5774, -0.5774, 0, 0, 0) and 
<p 2 = (-0.6969,0,0.6969,0.1196,-0.1196,0) « (-0.6969,0,0.6969,0,0,0). The two eigen- 
states are similar to the symmetry (ground state) and asymmetry (the first excited state) 
states of a single electron in three QD. This is also like the case in the driven two- level 
model. Moreover, the value of the parameters 2V/u = 2.4, which is the root of the zero- 
order Bessel function, suggesting that in this situation can be approximated by the driven 
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two-level model. 

We apply the perturbation approach^ to analyze the dynamical property for the strong 
field regime. When driven by strong field, the Hamiltonian(Eqinj) can be divided into two 
parts: H t which contain all the tunnelling terms, and Hi containing all interaction terms 
(these involving U\ and U2 and the electric field). We then find the eigen-system of the 
operator Hj(t) = Hj — i(d/dt), and employ the tunnelling Hamiltonian as the perturbation. 
When JJ\ = riiuj and U2 — ^1 the eigenvalue of Hi are 0. These represent the zeroth-order 
approximations to the quasi-energies in the perturbational expansion. The perturbation 
Hamiltonian H t induces the Floquet state transition. When the initial localization is |1), 
the nonzero transition matrix element is P x 5 ' = ((l|i?t|5 )) = s/2W J ni - n2 (V/u>). So the 
amplitude to frequency of driving field is the roots of Bessel function of n\ — ri2 order, the 
transition of initial localization state |1) is strongly depressed. By using progressive formula 
of Bessel function: J m {x) ~ cos(x — |m — |), we find the roots of J ni _ n2 are also the 
roos of \J2J^ 2 + 6Jni-n 2 when n\ — 112 is even. Therefore, at the point of quasi-energies 
collapse, the two-electron dynamical localization is set up. 

We show above the phenomena of dynamical localization occurs at the quasi-energy 
collapse when the two electrons are localized in left QD. This quantum phenomena can also 
happen when the two electrons are localized initially in centra QD. We verify this assertion 
by calculating the quasi-energy and BP min in Figure Efd) with the parameters U\ = 20, 
U 2 = 4, W = 1, u = 2, the initial condition C 2 (0) = 1, C a (0) = when a ^ 2. It shows 
that the features of the Figure Efd) are similar to that of the Figure Efb) for the strong field 
regime. The calculation indicates that the two electrons of initial localized in the centra QD 
can also maintain its localization at the quasi-energy collapse. For the weak field regime, 
the effects of Coulomb interaction emerge obviously. With the increase of U2, a series of 
quasi-energy crossing, and the dynamical localization is set up at the quasi-energy crossing. 

IV. SUMMARY 

In summary, applying three-site Hubbard model and Floquet theorem, we have investi- 
gated the dynamic behaviors of two electrons confined in a line-shape three quantum dot 
molecule driven by an ac-field. We find the triplet state and triplet state sub-space are 
decoupled fully, and can be investigated respectively. In special, the triplet state subspace 
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can be divided into 3 three-dimensional triplet state sub-space which are also decoupled. 
The analysis indicates the dynamical properties in the 3 three-dimensional triplet state sub- 
space, as well as the singlet state sub-space of no double-occupancy, are similar to that 
driven two electrons in two QD molecule. The calculations and analysis indicate that the 
dynamical properties in every subspace can be divided into two regime: In the strong-field 
regime, the dynamical behaviors are similar to that of non-interacting electrons in superlat- 
tice, and in the weak-field regime, the Coulomb interaction enhance the localization effect. 
We believe the method and the results presented here are useful to explore the effects of 
ac-field to control of multi-electron and multi-QD system. 
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TABLE I: The two electron distribution in three QD molecule 
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FIG. 1: The dependence of quasi-energies and (a) LP m i n , (b) CP m in as function of V, with the 
parameters: C/2 = 8, W = 1, u = 1. dotted line = quasi-energy exact results, dotted line = 
quasi-energy exact results, Solid line = perturbation theory. 
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FIG. 2: (a) magnified view of the weak field regime in Fig^a), (b) the quasi-energies and LP m i n 
as function of V, with the parameters: U2 = 16, W = 1, u = 1. dotted line = quasi-energy exact 
results, dotted line = quasi-energy exact results, Solid line = perturbation theory. 
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(a) (b) 




V V 

FIG. 3: The dependence of AP m i n as function of V with the parameters: (a) U\ = 0, U2 = 0, 
W = l,u = 2\ (b) Ui = 12, U 2 = 4, W = 1, u = 2. (c) magnified view of FigE^b) for the weak 
field regime, (d) The dependence of BP m i n as function of V with the parameters U\ = 20, U2 = 4, 
W = 1, u = 2. 



